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We study the ground state and excitation spectrum of a quasi one- dimensional nanostructure con- 
sisting of a pole and rungs oriented in the opposite directions ("centipede ladder", CL) at quarter 
filling. The spin and charge excitation spectra are found in the limits of small and large longitu- 
dinal hopping iy compared to the on-rung hopping rate t± and exchange coupling I±. At small 
t|| the system with ferromagnetic on-rung exchange demonstrates instability against dimerization. 
Coherent propagation of charge transfer excitons is possible in this limit. At large t\\ CL behaves 
like two-orbital Hubbard chain, but the gap opens in the charge excitation spectrum thus reducing 
the symmetry from SU(4) to SU(2). The spin excitations are always gapless and their dispersion 
changes from quadratic magnon-like for ferromagnetic on-rung exchange to linear spinon-like for 
antiferromagnetic on-rung exchange in weak longitudinal hopping limit. 

PACS numbers: 75.10.Pq, 71.10.Fd, 73.22.Gk 



INTRODUCTION 

The spin ladder systems attract much interest in both 
experimental and theoretical communities 1 during re- 
cent years. On the one hand, the ladder systems may 
be viewed as intermediate configuration between purely 
one-dimensional systems and periodic arrays of higher 
dimensions. On the other hand, in many cases the sim- 
ple theoretical models represent a good prototype for de- 
scription of fundamental phenomena like metal-insulator 
transition, formation of various density wave states, low 
dimensional superconductivity and many other strongly 
correlated effects 1 . Besides, some of the toy theoretical 
models are simple enough to be solved exactly on the lat- 
tice or in the long wave continuum limits. The impor- 
tant experimental realizations of the ladder systems (see 
e.gi 1 -) renewed an interest to imperfect spin chains and 
ladders affected by the chemical substitution, pressure or 
radiation defects. Although the literature on symmetric 
two-leg ladders at half-filling is quite extensive 1 -^, not 
so much is known about the effects of strong asymme- 
try in the exchange interaction (or hopping) between the 
legs^. Even less is known about the doped two-leg lad- 
ders away from the half-filling regime 6,7 . The most of 
considered cases are related to so-called incommensurate 
filling where in spite of strong Coulomb (Hubbard) in- 
teraction, the ground state of the ladder is metallic. On 
the contrary, at the commensurate filling (e.g. quarter 
filling) the Umklapp processes open a gap in charge ex- 
citations spectrum already for weak interaction. These 
effects can be consistently treated by means of bosoniza- 
tion technique^. 

In this paper we consider a model characterized by a 
gap opened due to kinematic constraints even before any 
interactions are taken into account. This model, which 
we call a centipede ladder (CL) or single-pole ladder (see 
Fig. la, c), can be visualized as the strongly asymmetric 



limit of the two- leg ladder (Fig. lb) where the interaction 
along one of the two legs is negligibly small compared to 
both interaction along the main leg and the rung, e.g. 
due to spatial orientation of the rungs (see Fig. [Tk.,b). 

The family of centipede ladder-type systems schemat- 
ically presented in Fig. QJ, was synthesized as a stable 
organic biradical crystal PNNNO^. The analysis of this 
geometry for the case of half-filling (two electrons on the 
rung with infinite on-site repulsion) has shown2ii£ that 
this model differs from the symmetric two-leg ladder. 
The half-filled case for the centipede ladder belongs to 
the specific class of universality manifesting itself in frac- 
tional (2/3) scaling of the spin gap. This new scaling was 
attributed to existence of hidden dynamical symmetries 
associated with quantum nature of spin-rotators formed 
at each rung of the ladder and should be contrasted with 
the integer scaling discovered for the symmetric two-leg 
ladders^. The corresponding model of highly asymmet- 
ric two-leg ladder where the interaction between the spins 
residing at the end of neighboring rungs is zero has been 
labeled as Spin Rotator Chain model^ because it may 
be mapped onto an effective ID chain with complicated 
on-site and inter-site exchange interaction revealing "hid- 
den" symmetries of spin rotator. 

One may expect that the same geometry of centipede 
ladder away from the half-filling may also demonstrate a 
behavior very different from such for symmetric two-leg 
ladder and/or the partially filled chains. A quarter-filled 
CL is the next commensurate system where these differ- 
ences should be looked for. At first sight, the quarter- 
filled CL, with strictly one electron per rung is equivalent 
to the quarter-filled two-orbital Hubbard chain (HC) 11 
with partially lifted orbital degeneracy. However, addi- 
tional local symmetries related to the spin and charge 
degrees of freedom of electrons on a rung influence a 
structure of the ground state and the spectrum of low- 
energy excitation and makes the phase diagram of CL 
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more reach than that of half-filled HC or quarter-filled 
two-orbital HC. The aim of this paper is to reveal these 
new features. The phase diagram of CL promises to be 
complicated enough, and in this paper we confine our- 
selves with several limiting cases. 

The manuscript is organized as follows: in Section I 
we formulate the model and describe the hierarchy of the 
electronic states on the rung. Section II is devoted to 
analysis of the limit, where the on- rung ferromagnetic 
exchange is dominant. In the Section III we consider the 
case of weak longitudinal hopping for antifcrromagnetic 
on-rung exchange. The Section IV is devoted to the case 
of strong longitudinal hopping and discussion of various 
response functions at quarter filling. Details of calcu- 
lations are presented in Appendix. In Conclusions the 
summary of results and the perspectives is discussed. 
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FIG. 1: (Color online) Centipede (a) and strongly asymmet- 
ric two- leg ladder (b). Several possible rung states (c) corre- 
sponding to vacancy, single occupied states, singlet and triplet 
states. Excited states with double on-site occupancy are not 
shown. 



I. MODEL AND HIERARCHY OF RUNG 
STATES 



The generic model Hamiltonian is H = H r + Hi , where 

H r = £ ^ Usl - a ~ tj - X! ( C L,o- C 0'> + H.c) (1) 
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describes the electronic states on the rungs, and 



Hi = 



V c t. 



C a i+l,<r + H.C. 



(2) 



is related to the longitudinal electron hopping along the 
leg. Here s = a, (3 enumerates sites belonging to the 
leg (a) and the branch in a given rung i. The 
Hamiltonian |2]) is is an extension of the Spin Rotator 
Chain Hamiltonian 9 ^ to the case where empty states on 
the rungs appear and electron hopping is possible both 
within a given rung (tj_) and between neighboring rungs 
(t||). We retained here the direct on-rung exchange in- 
teraction 7j_, which can be either ferromagnetic or anti- 
ferromagnetic. 

It is natural to use a single rung basis for characteri- 
zation of the energy spectrum of SRC. The set § of 16 
rung basis states consists of an empty state |0) (which 
we will call a vacancy), four singly occupied rung states, 
four doubly occupied rung states which are split into sin- 
glet/triplet configurations, two doubly occupied on-site 
states, four triply occupied states and a fully occupied 
state with four electrons. In the limit of strong Hubbard 
repulsion, U 3> , the lowest states on the rung i belong- 
ing to the charge sector iVj = 2 (doubly occupied rungs) 



|2-5,i) = -^(4 T ct. i 



|2T1,«) 
|2T0,t) 
\2Tl,i) 
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(3) 



Here the quantum numbers of two-electron states are 
qualified as a singlet 25 and triplet 2TA with A = 1, 0, T. 

The energies corresponding to the eigenstates (J3|) arc 
defined as follows 



Ei,u — e o, E 2 s,i — 2e — j, EiTi — 2eo — I A 



(4) 



where | cq I is a single-electron ionization energy, j = 
2t 2 L /U is the indirect exchange of spins at adjacent sites, 
which stems from virtual transitions including double oc- 
cupancy. Strong Hubbard repulsion effectively excludes 
all states beyond the four states ^ in the low-energy 
part of the excitation spectrum with Ni>2 (see below) . 

In a charge sector Ni = 1 the electron spectrum de- 
mands further diagonalization. The states on a given 
rung are characterized by spin and pseudospin. Spin Sj 
is defined as usual, and pseudospin n is introduced as 



2 ^2 [ c /3^ c 



H.c 



(5) 

The kinetic energy term in ITJ corresponding to the hop- 
ping along the rung may be cast in the form 



K, I = — 2t\ 



(6) 



while, assuming two different electro-chemical potentials 
fi a and [10 (or, equivalently, two different surroundings 
for a,P sites), we get one more term 

^Tifaa-llti-^^NiQic+Hf,) (7) 
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in H r - Here Ni — 7ii a +riip. Terms in ©,((21) proportional 
to t x and t z correspond to external uniform " magnetic" 
field, causing the Zeeman splitting of the rung orbital 
states. We discuss the formation of the orbital gap in 
the Section IV, taking everywhere below /i Q = (ip. 

If the sites a and (3 are equivalent, the non-diagonal 
transitions described by rf are eliminated in the local 
basis of the bonding (+) and anti-bonding (— ) states 
introduced in accordance with 



\±,a,i) 



1 



(8) 



At finite tj_ this basis still diagonalizes the rung Hamilto- 
nian and the corresponding single-electron energy levels 
are e± = e T t±. One should stress that the doubly oc- 
cupied states cannot be viewed as two bonding and/or 
antibonding orbitals on a rung, due to strong Hubbard 
repulsion, and the charge sector Ni — 2 is represented 
by the states ([3]). This is a specific property of CL in 
comparison with the two-orbital HC. 

Another remarkable property of CL which, in fact, al- 
lows double occupation of a rung without paying large 
energy U, is its instability against chain dimerization un- 
der certain conditions. To describe this instability, one 
should introduce the elementary cell for a CL with dou- 
bled lattice period. Such cell is occupied by two electrons, 
and its energy spectrum at zero tii is easily found. The 
low-energy subset is formed by the dimerized states and 
the states with singly occupied rungs, respectively: 

{Ni = 2,N i±1 = 0} = {|2rA,t;0,t + l),|25,»; 0,*+l)} 
{Ni = l,N i+1 = l} = {\±,a,i;±,a',i + l)}. (9) 



The corresponding two-electron energy levels are 

£r(2,0) = 2e -/x, 
£J<j(2,0) = 2eo- j, 
£ ++ (l,l) = 2(e -tx) 
B+_(l,l) = 2eo 
B__(l,l) = 2(e +t±) 



(10) 



One concludes from (|T0|) that the lowest energy levels 
which predetermine the spectrum of quarter-filled CL are 
E T (2,0), E s (2,0), £++(1,1), ^+-(1,1). If I± < 2t±, 
the lowest two-electron states of this elementary cell be- 
long to the subset Sh = {+,<r,i; +,a' ,i + 1} with one 
bonding orbital per rung. In the opposite case I± > 2t± 
the lowest state correspond to (2,0) occupation of ad- 
jacent rungs. We denote the corresponding subset as 
§d = {2T, 2S, i; 0, i + 1}. At finite longitudinal hopping 
CL forms homogeneous quarter-filled chain in the former 
case and dimerized chain with alternating (2,0) occupa- 
tion in the latter case. 

In case of negative (antifcrromagnetic) exchange I± 
the ground state is homogeneous quarter-filled chain be- 
longing to the subset §>h — {+, a, i; +, a ', i + 1} and the 
lowest excited state is E$ (see Fig. Wp)- 
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FIG. 2: (Color online) Two electron level structure for two 
neighboring rungs, a) Ferromagnetic on-rung exchange. The 
ground state is dimerized. b) Antifcrromagnetic on-rung ex- 
change. The ground state is array of bonding orbitals. c) Oc- 
cupation scheme for two different ground states: filled, empty 
and half-filled circles correspond to occupied (triplet), empty 
(vacancy) and bonding electron states respectively. Transi- 
tions within a given charge spin or charge sector and transi- 
tions between different sectors are marked by vertical and 
slanting dashed arrows, respectively. The energies of this 
transitions are also shown. 



In two following sections we discuss charge and spin 
excitation spectra of CL in the limit of weak longitudi- 
nal hopping. Looking at Fig. [^,b and Eqs. (fTT)]). we 
see that in case (a) one may find a regime, where only 
the levels E ++ and Et are involved in formation of the 
low-energy excitations, whereas in case (b) at least three 
states £++, E$, E^ are involved. 



II. WEAK LONGITUDINAL HOPPING LIMIT, 
FERROMAGNETIC ON-RUNG EXCHANGE 

In this section we analyze the quarter-filled CL with 
7+ > in a situation , where the longitudinal hopping is 
small in comparison with the on-rung coupling constants, 
and the inequalities 



I\ >2t 



*H < 1+ - 2t A 



(11) 



are valid. In the absence of longitudinal hopping each 
second rung is empty and all doubly occupied rungs are 
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in spin one triplet states (see Fig. [2^). One has to 
study longitudinal spin and charge excitations at finite 
tii above the dimerized ground state with alternating 
rungs . . . |2TA, i), |0, i + 1) . . . The charge excitations arise 
due to transitions E>d Sft. These excitations (charge- 
transfer excitons) are obviously gapped, so that the low- 
energy excitations belong to the spin sector. 

In zeroth order approximation the spectrum of rung 
dimers is determined by the eigenstates of the Hamil- 
tonian fl} in a charge sector Ni = 2 for FM on- 
rung exchange. These dimers are in the state of spin 
rotator ^ 10 ]13 characterized by the singlet-triplet mani- 
fold {|2TA,i), \2S, i}}, with the energy gap between two 
spin states given by 



AE 



ST 



I\ 



(12) 



Longitudinal hopping ^ generates indirect exchange 
produced by the electron cotunneling via empty odd 
rungs. To derive the corresponding spin Hamiltonian, 
we have to exclude the charge sector Sh from the ef- 
fective phase space. It is seen from Fig. [2^,, that 
the singlet states are effectively quenched at the energy 
scale I± — 2t±, which characterizes the charge transfer 
gap, so that one may confine oneself with the subspace 
{S^jS/,.}, with singlet states excluded from the subset 

= {2T, i; 0, i + 1} of "polar" states. 

The remarkable property of the state E>' d is that the 
macroscopic configurational degeneracy of the Hubbard 
model with partial filling is removed due to dimeriza- 
tion. The only remaining even-odd degeneracy of the 
dimerized ground state of zero-order Hamiltonian may 
be also removed in case of CL with odd number of rungs 
M = 2N + 1, where, say, all even rungs 0, 2, . . . 2N are 
doubly occupied and all odd rungs l,3...2iV — 1 are 
empty. Then only the trivial spin degeneracy remains, 
and one may strictly define the effective exchange Hamil- 
tonian in the 4-th order of Brillouin-Wigner expansion 



Keff = % 



* |h, (4) i* 



(13) 



where the brackets (^q\ • ■ ■ l^o) stand for the dimerized 
polar ground state with alternating triplet spin states 
and vacancy states 



(4) _ H^ (1) )(0 (1) |^|0 (2) )(^ (2) I^I0 (3) )(0 (3) |^ 



Hl (-4 X) )(-4 2) )(-4 3) ) 

Here the projection operators 

p(oQ = !«/,(«)) |, a = 1,2, 3 



(14) 



include the states which are generated from the ground 
state l^o) after 1-st, 2-nd and 3-rd action of the operator 

Hi, respectively and Eq** are the corresponding excita- 
tion energies. The 4-th order cotunneling process result- 
ing in effective kinematic exchange is illustrated by Fig. 
02 It is evident that only the clusters (i, i±2) are in- 
volved in formation of the effective exchange Hamiltonian 




FIG. 3: Scheme of superexchange induced in 4-th order hop- 
ping processes in dimerized CL. Arrows indicate sequential 
electron hops resulting in the intermediate states (|15[) . No- 
tations for the occupation scheme are the same as for Fig. 
2c. 



where i stands for any even site. We consider here the 
simplest situation, where the singlet intermediate states 
can be discarded in the main order, which happens at 
t\\/(I± — j) <C f||/(ij_ — 2tjJ. This simplification does 
not change our qualitative conclusions. 

= |+o-i,i; + a 2 , i + 2T\,i + 2) 
\^) = \+a u v t 2TA 1) i+l; +a 3 ,i + 2) (15) 
|0 (3) ) = |+<7i,i; +CT 4 ,i + l; 2TA',i + 2) 

It is seen from Eq. (fT5|) and Fig. [3J that the intermedi- 
ate states involve electron cotunneling exchange between 
two nearest even sites with formation of a virtual dou- 
bly occupied triplet state on the odd site between them. 
This mechanism is known as superexchange in the the- 
ory of magnetic dielectrics. Appearance of intermediate 
triplet states |2TA) in the indirect exchange interaction 
is the peculiarity of CL: unlike the two-orbital HC model, 
double occupation of intermediate state does not cost the 
Hubbard repulsion energy U , because the two electrons 
on a rung may occupy different sites. As a result we come 
to the following effective spin Hamiltonian 



K 



eff 



J||J> 



(16) 



where the site index i stands only for even rungs of CL. 
The coupling constant is 



'11 = 



1 1 



16A| 



(17) 



where Ax = ij_ — 2t± is the charge transfer energy gap 
(see Fig. [5^). The details of calculation are provided 
in Appendix. Deriving (|17[) we assumed that the sin- 
glet state does not contribute to superexchange. This 
assumption holds provided A^ «C AE$t- 

Thus we have demonstrated that the problem of 
quarter-filled CL is mapped onto the familiar problem of 
spin-one Heisenberg chain with ferromagnetic exchange 
in the limit of weak longitudinal hopping under condi- 
tions (|11|) . Only the even sites are involved in formation 
of the gapless spin excitations, which, apparently may be 
described in terms of the spin-wave theoryJ^ 
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FIG. 4: (Color online) Tight-binding band structure for quar- 
ter filled centipede model at tu <C t±. The occupied and 
empty orbital subbands e±(fc) are shown. The chemical po- 
tential is chosen as a reference energy level. 



As to the charge transfer excitons, | + a, i; +a' , i + 1), 
shown in the right column of Fig. [2J;, these states, like 
other excitons can exist both in triplet and singlet spin 
configurations. These excitations may propagate coher- 
ently by means of cotunncling reminding that described 
byEqs. ITU) . (fH))) . To provide translation of such exciton 
from one cell of dimerized CL to another, two electrons 
should hop to the right (i —> i + l;i + 1 —> i + 2) and 
then to the left (i + 2 -> i + 1; i + 3 -> i + 2). This 4-th 
order indirect "exchange" provides excitons with the dis- 
persion law uj q ~ Ay + 2Ju cos 2qa, where a is the lattice 
spacing along the leg. 

Thus, we have shown in this section that the spin ex- 
citations in the CL with ferromagnetic on-rung exchange 
are gapless magnon-like modes, whereas the excitations 
in a charge-sector are gapwise excitons describing coher- 
ent propagation of two bonding states along the chain. 



III. WEAK LONGITUDINAL HOPPING LIMIT, 
ANTIFERROMAGNETIC ON-RUNG 
EXCHANGE 

In case of negative on-rung exchange hopping along the 
leg and small ty, so that the inequality 



t|| < 2t ± - j 



(18) 



is satisfied, the energy levels are ordered in accordance 
with Fig. r^b, and the ground state wave function of 
quarter-filled CL is given by the product of bonding or- 
bitals 



l*o> 



1; 



(19) 



In this sector of phase diagram, the quarter-filled CL is 
close in a sense to the quarter-filled two-orbital HL, 11 
with lifted orbital degeneracy: each rung is occupied by 
one electron with bonding orbital, and the empty anti- 
bonding level is separated from this band by the orbital 
gap - 2t±_ [see Eqs. ©, © and Fig. 0]. 

The important distinction from two-orbital HL roots 
in the structure of the doubly occupied states. Since two 
electrons on a rung occupy different orbitals in accor- 
dance with the energy cost of this state is not the 
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FIG. 5: Coherent motion of the vacancy-singlet exciton due to 
the 4-th order electron hopping (1 — » 2; 2-^3; 3 — > 2; 4 — > 
3). Diagrams illustrate electron configurations in the initial 
(a), subsequent intermediate (b) — (d) and finite (e) states 
of exciton translation. Filled, empty, half-filled (black down) 
and half-filled (black up) circles correspond to singlet, vacancy 
bonding and antibonding electron states, respectively. 



Hubbard repulsion energy U but the difference As = 
Es — £■++ = 2tj_ — j, which is even smaller than the 
orbital gap. Nevertheless, magnetic properties of these 
two models are similar in the limiting case (fT5)) : the lon- 
gitudinal hopping generates effective AFM longitudinal 
exchange ~ ijj/A,g, and the gapless spin liquid state of 
the resonance valence bond (RVB) type is realized. 

The charge excitations exhibit more interesting char- 
acter, as we discuss below. The lowest excited state is a 
vacancy-singlet pair occupying neighboring rungs 

\S,i;0,i+l) = \ ... + a,i-l; S,i; 0,i+l\ +a,i + 2...) 

(20) 

with the energy Ag. 

Clearly, the state 15, i; 0, i+1) is an analog of the charge 
transfer exciton briefly discussed in the end of section II. 
However, in this case the exciton propagates via inter- 
mediate bonding-antibonding levels e± due to the strong 
inequality t± ^> j. Figure [5] illustrates the mechanism of 
exciton propagation. Here the rung occupied by the anti- 
bonding electron is shown by the pair of half-filled circles 
with filled upper half. Then the elementary hopping from 
doubly occupied to empty rung may be presented as 



\S,i;0,i + l) -> \±,i;T,i + l) 



(21) 



and it costs the energy j, according to (fTO]) . Note that a 
crucial detail allowing one to discuss exciton propagation 
as the charge excitation is the fact, that the spin sequence 
of the bonding orbitals in the initial and final state in Fig. 
[5] is unchanged and the spinon dispersion is unaffected 
by the coherent motion of the singlet exciton in the long 
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wave limit. This is a manifestation of the spin charge 
separation in quarter-filled CL model. 

In the limit of extremely weak longitudinal hopping 
t\\ j, the singlet-vacancy pair may propagate coher- 
ently. This pair is doubly degenerate due to its permuta- 
tion symmetry, and this degeneracy manifests itself in ex- 
istence of right- and left-moving charge-transfer excitons. 
The 4-th order hopping process resulting in translation 
of the right-mover between two adjacent double cells is 
shown in Fig. [5] Superposition of right- and left-moving 
excitons results in the dispersion law 



As — J'\\ cos 2qa 



(22) 



with Jjj ~ t^/j 3 - We showed before that the spin se- 
quence remains unchanged (Cf. RefH3) upon the motion 
of exciton, and the gain in the kinetic energy is ~ Jjj . The 
resulting spin configuration in Fig. [5] corresponds to two 
spins detached from the main spin sequence on the right 
and attached to one on the left. For long-wave modes of 
exciton motion, q — > in (|22|) , the change in the mag- 
netic energy associated with this disruption is negligible 
compared to the gain ~ Jjj , and charge excitons move 
coherently. 

Now we turn to the case t» > j. In this situation the 
coherent "vacancy-singlet" excitons are unstable against 
dissociation into independently propagating singlet and 
vacancy " defects" due to the overlap of the bands orig- 
inating from the levels E$ and Apparently, only 

incoherent propagation of such defects is possible. An ex- 
ample of 4-th order process resulting in decay of singlet- 
vacancy pair is presented in Fig. [5J Such incoherent 
processes will be accompanied by spin-flips, antibonding- 
bonding de-excitation and eventual relaxation towards 
the ground state. 

However, coherent motion of defects is possible at small 
deviation from one-quarter filling. Excess holes (vacan- 
cies) behave similarly to the holes in the nearly half-filled 
HC in the limit (JTHJ) , where the longitudinal hopping does 
not result in excitation of the upper orbital subband e_. 
In this case one deals with a fermion gas, "fully polarized" 
with respect to the pseudospin r. Excess electrons fill the 
side sites on the rungs without paying the Hubbard en- 
ergy U, and it is easy to check that the hopping along 
the leg is not accompanied by creation of a spin string 
because the doubly occupied state is a spin singlet. As 
a result, the coherent propagation of these defects is also 
possible under the conditions (fT5|) , when the antibonding 
orbitals are not excited. Due to specific geometry of CL, 
there is no particle- hole symmetry near the 1/4 filling, 
and the hopping rates for particle and hole propagation 
are different. Direct calculation shows that the hopping 
rate for coherent motion of a vacancy reduces due to the 
kinematic constraint from fn to t v — £ii/2, whereas the 
similar reduction for the singlets gives t s — t^/A. Thus, 
basically the coherent motion of the vacancy and sin- 
glet can be considered in terms of tight-binding bands 
£ i,2(fc) = — 2ti ]2 cosfc. Evidently, singlet (vacancy) can- 
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FIG. 6: Incoherent motion of the vacancy and singlet due to 
the 4-th order electron hopping (2 — » 3; 2 — > 3; 1 — > 2; 3 — > 4). 
Notations for the occupation scheme are the same as for Fig. 5. 



not occupy the same rung twice, hence the corresponding 
excitations are fermion-like. We can think of them as of 
spinless fermions characterized by two colors fi^ v and fi. s : 



|0,i><+,*|-/L \S,i)(+,i\^fl 



(23) 



Even at exact quarter filling but at finite temperatures, 
the vacancies and singlets may appear at small concen- 
trations due to thermal activation. These excitations are 
described by the effective two color Fermi liquid model 



H, 



eff 



(24) 

where the effective constant K describes effective short- 
range repulsion between two fermions, which prevents 
them from occupying the same site. This constant may 
be estimated as K ~ tjj/Ag, which is small comparing 
to the band- widths ~ in. Having in mind the asymmet- 
ric hopping Hubbard model 20 , it is tempting to analyze 
it using the well-established machinery of bosonization 
etc. However, the sizeable concentrations of vacancies 
and singlets arc expected at temperatures comparable to 
the main parameter, t±. This argument shows that that 
fermion gas in H e tt is not degenerate, and the lineariza- 
tion of the spectrum around Fermi points is unjustified. 

As to the excitations in the spin sector, one may say 
that the problem of spin excitations in presence of a sin- 
glet or a vacancy may be mapped on the problem of RVB 
spin liquid with spin-one defect or dangling bond, respec- 
tively (see Ref. [lg for general approach). 

To summarize, we have found in this section, that the 
spin excitations are gapless Fermi-like spinons both at 
exact 1/4 filling and at small deviation from this point. 
As to the coherent charge excitations, these are gapped 



excitons at extremely small tn < j and gapless fermions 
at finite doping and/or finite temperature. 



IV. STRONG LONGITUDINAL HOPPING 

The CL model possesses unusual properties also in the 
limit of strong longitudinal hopping, when it dominates 
over the transversal degrees of freedom: 



ill > tx » Ii 



(25) 



In this case the on-rung excitations shown in Fig. [2] 
are completely smeared, and the structure of low-energy 
spectrum of the system is predetermined by orbital de- 
grees of freedom (bonding-antibonding subbands) and 
corresponding Hubbard repulsion. We will show in this 
section that in this simple situation close to the two- 
orbital HC the charge excitations behave unconvention- 
ally because of the gap opening and corresponding van 
Hove singularities appearing in the lower partially occu- 
pied Hubbard band. 

For a clear presentation we compare the case of the 
asymmetric two-leg (2L) ladder and centipede ladder (see 
Fig. [lb, a, respectively. Fourier transforming along the 
leg, we write the kinetic part of the Hamiltonian as 



Hki 



2t\\ cos fc 



■fi, t± 

2Pt\\ cosfc — ft 



*k ~ \ c a,k,a' f -B.k,a) 



*fc(26) 

(27) 



where P = 1 for the 2L model and P = for CL model. 
At this point we also allow for the difference in the elec- 
trochemical potentials, 2/i = fi a —fig, between the sites a 
and ft [see Eq. (JJJ)]. The dispersion in these two extreme 
cases is given by 



e±(k) = —2i|i cosfc q= y/J, 2 + t 



2L 



(28) 



= -t\\ cosfc T y*i + (A* + *|| cos A;) 2 , CL (29) 

The square root term in Eq. [28] has a meaning of effec- 
tive Zeeman splitting. In contrast to 2L situation where 
the bands e_|_(fc) and e_(fc) bands overlap we see that in 
CL case two branches of dispersion are always separated 
by the gap (Fig. [7]). This indirect gap is given by 



A ind = 2(^*2 +i| -t||) ~t?_/t||. ( 30 ) 
for the case of equal electrochemical potentials (/i = 0) 



A ind ^ti*||/(tjj -fi 2 ) 



(31) 



for fi a 7^ fig. The direct gap is realized at fc = 7r/2 
for fi = and at fc = arccos(— fi/t\\) otherwise and is 
equal to A = 2t± in both cases. For exact 1/4 filling, 




FIG. 7: (Color online) Tight-binding band structure of two 
lowest Hubbard subbands for quarter filled centipede model 
at t± <C t||. The electro-chemical potential fi a = fib- The 
direct gap corresponds old Fermi surface point &f = n/2. 
The indirect gap Ai n d is between the bottom of conduction 
band and the top of the valence band is characterized by the 
"flat" parts of the spectra. The inset shows the energy de- 
pendence of the Density of States (DoS). The character of the 
DoS singularities and its influence on the response functions 
is discussed in the Section IV. 



the overall chemical potential fi$ = fi a + fig lies within a 
gap separating two bands. The existence of two (direct 
and indirect) gaps makes the CL case principally different 
from 2L model. 

One could expect that in the limit tii S> ij_ the proper- 
ties of the system will be close to those of two-orbital 
Hubbard chai n 11 ! 21 because the orbital degeneracy is 
nearly removed by strong longitudinal hopping. This is 
not the case, however. The gap in the spectrum prevents 
formation of ST/ (4) manifold formed by spin and orbital 
degrees of freedom. However the interband transitions 
as well as the edge singularities in the DoS influence the 
thermodynamic and optical properties of CL 

In the remainder of this section we concentrate on 
the case fi — and consider the manifestation of sin- 
gularities in the excitation spectrum in different re- 
sponse functions of CL. The density of states p{E) — 
(27r) _1 j dkS(E — e(fc)) is divergent at E = ±t±_ and at 

±(*l| + Jt 2 , +t 2 ). It behaves as 



E 



P(E) 



A 



-1/2 X^$( 
d 



-Ajnd ± E) 

V ±E — A in d 



(32) 



at small energies \E\ 



here = 1 at x > 0. 



Such character of the density of states (see inset in Fig. 
[7]) should lead to observable anomalies in the thermody- 
namic quantities. The specific heat C, compressibility 
X and dc conductivity cr^c are sensitive to indirect gap. 
Setting all fundamental constants to unity we write for 
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thermodynamic quantities per unit cell 
E 2 dE 



C 



4T 2 coslr (£/2T) 



P(E) 



T 
T 



1/2 



T> A 



ind 



3/2 



e -A Ind /T T <<; ^ (33) 



-Am 



for the specific heat and 

f dE 

X ~ J 4Tcosh 2 (£/2T)' 

~ (TA md y 1/2 , T>A ind 

~ (A itK ,T)- 1/ V A '»« , / r , T«A md (34) 

for the compressibility along the ladder. Assuming the 
elastic quasiparticle lifetime r e ;, the dc conductivity is 
given by 



<?dc = T e l 



dk f de ± (k) Y 

4Tcosh 2 (e ± (fc)/2T) \ dk J 



~ T el ^TA md , T»A 

~ Tei V / TA ind e- Ai "< i / T , T<A md (35) 

In contrast to behavior of thermodynamic quantities 
sensitive to indirect gap, the optical conductivity re- 
sponse is sensitive to direct gap. It is non-zero at \lo\ > A 
and is given by^ 2 -: 



Copt (ui) 



- I tflfe[n(e_(*)) - n(e+(k))} 



t\t 2 sin 2 k 
t\ + t 2 cos 2 k 
i|,A 2 



u} 2 Vw 2 - A 



i _ A 2 
iiiA 2 



S(u + e_(fc) - e + (fc)) 
T< A 
. T > A 



(36) 



We thus see a continuous absorption band in cr opt (o;) in 
contrast to the case of 2L ladder, where the optical con- 
ductivity consists of a single line, a op t{oj) oc S(oj ± A). 

Interband optical transitions may result in creation of 
excitons with a center-of-mass momentum Q = %. Un- 
like the case of weak longitudinal hopping (section II), 
these excitons arise in the conventional way due to attrac- 
tive Coulomb interaction between electrons and holes. 
The bound electron-hole states formed on the almost flat 
parts of the spectrum over the indirect gap are character- 
ized by heavy effective mass. We therefore expect that 
these excitons are nearly localized and one may neglect 
the processes of their coherent motion. 



V. CONCLUSIONS AND PERSPECTIVES 

In this paper we have demonstrated several fragments 
of the rich phase diagram for centipede ladder. It was 
shown in our previous studies2J£ that the half-filled CL 
is an example of spin 1 chain with soft triplet-singlet 
excitation, which does not belong to the Haldane gap 
universality class. Now we have found that the 1/4 filled 
CL also demonstrate unconventional properties both in 
spin and charge sectors. 

We considered here only three limiting cases of the 
quarter filled CL where the character of charge and spin 
excitations may be revealed by means of relatively simple 
arguments. First, we have found the regime where CL 
behaves as a spin 1 chain with ferromagnetic effective 
longitudinal exchange and magnon-like excitation spec- 
trum (Section II). It is shown that the system is unstable 
against dimerization at strong enough on-rung ferromag- 
netic exchange. The only coherent mode in the charge 
sector is the gapped triplet charge-transfer exciton. 

One should stress that the dimerisation effect offered 
in this section has no analogs in current literature. Usu- 
ally the dimerization mechanism in spin chains is ei- 
ther explicitly built in the Heisenberg Hamiltonian or 
arises as a result of competition between the nearest and 
next nearest exchange coupling 2 ^. As to the two-leg lad- 
ders, the dimerized configurations were considered in a 
context of RVB states (rung-dimer states vs leg-dimer 
states In our case the lattice period is doubled as a 
result of competition between the on-rung hopping and 
exchange coupling. Experimentally, rung dimerization as 
a result of charge redistribution between legs and rungs 
was detected optically^ in two-leg compounds AVgOi5 
(A=Na,Sr). However, in vanadium bronzes this effect is 
driven mainly by the valence instability of V ions. 

Second, we have described the situation, where the va- 
cancies (empty rungs) and spin singlets (doubly occu- 
pied rungs) may propagate either as coherent excitons or 
as independent spinless fermionic quasiparticles with dif- 
ferent Fermi velocities (Section III). The dispersion law 
is gapped in the former case and gapless in the second 
case, where the charge sector may be described as a two- 
component Fermi liquid with weak short-range interac- 
tion. The spin subsystem behaves as the RVB-type spin 
liquid, where the spinons propagate in presence of spin 1 
and dangling bond defects. 

Third, we have considered the limit of strong longitu- 
dinal hopping, where the CL behaves as a specific version 
of two-orbital Hubbard chain. The peculiar features of 
this limit is the appearance of a gap in the half-occupied 
lower Hubbard band. This gap has purely hybridization 
nature, but the van-Hove singularities in the density of 
states strongly influence the thermodynamic and optical 
properties of CL. 

Our analysis of these three cases by no means ex- 
hausts the variety of unusual properties of CL at partial 
filling. In particular, we didn't include in our Hamil- 
tonian direct longitudinal exchange which can enhance 
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the indirect kinetic exchange or compete with it. We 
didn't consider the phase transitions between different 
phases, which may have their own peculiarities as com- 
pared to currently studied quantum phase transitions in 
low-dimensional systems 2 ^. We also didn't discuss pos- 
sible lattice distortions, which may accompany dimeriza- 
tion. All these and many other open questions are left 
for future investigation. 

The centipede ladder shown in Fig. [1] is the simplest 
example of the family of quasi ID systems intermedi- 
ate between chains and two-leg ladders. The experimen- 
tal realization of such system (organic biradical crystals 
PNNNO 8 ) is already known, but one may easily imagine 
molecular chains decorated not with rungs but with rad- 
icals forming closed loops, zigzags etc. We believe that 
the study of excitation spectra, magnetic and transport 
properties of this family will bring new unexpected re- 
sults. 
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which leads to Hf$o,o = — "4$o,o- Similarly, one obtains 

t 2 ~ 

$2,2 = -9^2, 2- Because the projected action of Hf is 
self-conjugate, we obtain Hf $0,0 = 16*0,0 an d Hf $2,2 = 

4*2,2. 

We would also like to exclude more complex structure 
of the effective Hamiltonian, namely 

H *ff = - J \\ E( S * ' S *+ 2 + Cl ( S * ' S *+ 2 ) 2 ) ( A - 2 ) 

i 

which is not prohibited by symmetry. 

To this end, we calculate also the correction to states 



APPENDIX: CALCULATION OF THE 
SUPEREXCHANGE INTERACTION 

In order to calculate J11 we note that two spins S = 1 at 
even sites in Eq. (|16p and in Fig. [3] can be combined into 
total spin S to t = 0, 1, 2. The scalar product Si • S^+2 = 
^((Sj + Sj+s^-Sf-S 2 ^) in (16) attains values -2, -1, 1 
for these values S to ti respectively. The value of 3J|| can 
hence be naively calculated as a difference in fourth-order 
corrections (TT4")) to the states Stot = and Stot = 2. 

Let us use a shorthand notation |m, m!) for a quan- 
tum state \Tm,i;Tm',i + 2), and <i>L,m for the state 
with total spin Stot — L and its z— component m. The 
weight of $L,m" in |m,m') is the Clebsch-Gordan co- 
efficient Cf r J™ lm , . Then the states we are interested in, 
Stot = and S to t = 2, are given by $0,0 = (| - 1, 1) + 
|l,-l)-|0,0))/V3and $ 2 , 2 |1,1). 

Further, we denote the combinations i\ +03, i+2) 
in |</>( 2 )), Eq. (fl"5)) as "non-local triplet" states, e.g. 

|+ T,i;+ t,« + 2} = |?M + 1) etc. The interme- 
diate state in Fig. 3, obtained after two longitudinal 
hops in terms of its spin components, is rather obvi- 
ously classified in the same terms of total spin, com- 
posed of two triplets, \Tm, i + 1) and \Tm, i + 1). Simi- 
larly denoting \Tm, i + 1; Tttt/, i + 1) = \m, ml) we have 
$0,0 = (hU 1) + 1 1,-1) - |0,0))A/3 etc. 

It then can be shown that (projecting at each step at 
the states with lowest energy, bonding and triplet, as 
discussed above) 



(A.l) 

I 

$i. m of total spin L = 1. Within the same logic, we have 

wf(|i,-i)-|-i,i)) = HfV2$i 1 o - f |s,; + i ; T(M+i) 

where S non-local singlet. Therefore Hf<&\fl = |i|$i,o 
and for our values L = 0,1,2 we have the corrections 
-1/16,-1/8,-1/4 in units of ijj/A|, respectively. It 

means that J|| = tjj/(16A|) and c x = 0. Thus the 

anisotropy term (S-; • S.; + 2) 2 does not appear in the forth 
order of longitudinal hopping. The reason of it is very 
simple: in order to get the four-spin S = 1 interac- 
tion within the scheme shown on Fig. [3] where four spins 
s = 1/2 participate, one needs at least 8 processes of 
longitudinal hopping. 



Wf|-l,l) - J=|+|,i;T0,i+l;+T,i + 2) 
Hf(\ - 1, 1) + |1, -1)) = ||070> = |[yf $2,0 - y^o.o] 
«?|o, o) = i> + |i, -1) + |o, o)) = -$ 2 ,o + -$ ,o] 
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